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ALGORITHM FOR FUEL CONSERVATIVE HORIZONTAL CAPTURE 

TRAJECTORIES 

Frank Neuman 
and 

Heinz Erzberger 
SUMMARY 

This paper describes a real-time algorithm for computing 
constant altitude fuel-conservative approach trajectories for 
aircraft. It is a refinement in the horizontal plane of thw 
fuel conservative capture algorithm previously flight tested 
at the Ames Research Center. The trajectories may start at 
any initial position, heading and speed and end at any other 
final position heading and speed. They are geometrically 
simple, as they consist only ot straight lines and a series of 
circular arcs of varying radius. Throttle control is also 
simple as it consists of maximum thrust in the acceleration 
segments of flight, nominal thrust to balance drag in the 
constant speed segments, and zero thrust in the deceleration 
segments. The bank -angle control is either zero or ' 
approximately 30°. The algorithm is designed to approximate 
constant bank -angle decelerating turns by combining a series 
of short circular arcs of different radii into a single turn. 
The bank angle varies slightly due to flying the circular arc 
segments with varying speed. 


To study possible further fuel reduction, a more sophisticated 
algorithm, which is suited for long-range flightpaths but not 
short-range paths, was developed that uses constant bank-angle 
turns (rather than a sequence of circular arcs to approximate 
a constant bank angle) . The more sophisticated algorithm is 
too slow for onboard implementation and was only investigated 
to serve as a basis of comparison with the real-time 
algorithm. 

The characteristics of the trajectory computed by the two 
algorithms were chosen to approximate the extremal 
trajectories obtained from the optimal control solution to the 
problem. The fuel consumption of the trajectories computed 
with the two algorithms was compared with the corresponding 
of the extremal. For the most frequently occuring 
trajectories in the terminal area the comparison with the 
extremals showed a fuel difference of only 0.5 - 2% for the 
real-time algorithm in favor of the extremals. Zero- thrust 
extremals were the most difficult to approximate by the 
algorithm and the synthesized trajectories took up to several 
hundred pounds more fuel than the corresponding extremal. 
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INTRODUCTION 


The escalating costs of fuel have increasingly focused 
attention on the problem of minimizing fuel consumption of 
commercial aircraft. Recent work in aircraft guidance has 
demonstrated that onboard optimization of aircraft 
trajectories offers an efficient method for fuel conservation. 

Aircraft trajectory optimization problems can be divided 
into two classes, namely enroute problems with trajectory 
lengths of 50 n.m. and longer, and terminal area problems with 
lengths of from three to 50 n.m. For the enroute problem, an 
onboard algorithm for fuel/cost optimum climb-cruise-descent 
has been developed by application of optimal control theory 
(ref. 1). The solution was developed for the vertical plane 
only, since horizontal maneuvers are not a significant feature 
of enroute flight. 

Terminal area trajectory problems are more difficult to 
solve since vertical and horizontal maneuvers involving speed, 
altitude and heading changes occur simultaneously, and are of 
comparable significance in influencing fuel consumption. The 
implementation of the exact solution to the optimal control 
problem is too complicated for real-time mechanization in a 
flight computer. As a result it is necessary to develop 
approximate solutions to the problem. An algorithm for 
onboard calculation of fuel conservative terminal area 
trajectories was recently described by Erzberger and McLean 
^ef. (2)). Their algorithm computes fuel conservative approach 
trajectories, which they refer to as capture trajectories, 
starting from an arbitrary initial position, heading airspeed 
and altitude and ending at a specified final position, 
heading, airspeed and altitude. The final position is chosen 
as close to the touchdown point as possible. 

The algorithm of Erzberger and McLean separates the 
problem into two parts in order to achieve computational 
simplicity: first, the synthesis of a horizontal trajectory 
consisting of constant radius turns (up to three) and a 
straight line; ana second, the synthesis of fuel conservative 
airspeed and altitude profiles along the previously computed 
horizontal trajectory. The synthesized horizontal trajectory 
minimizes the distance to fly. 

It is known from several studies of fuel and time optimum 
trajectories (ref. (3), (4), and (5)) that horizontal maneuvers 
flown during simultaneous accelerations or decelerations often 
require turns at maximum bank angle rather than turns at 
constant radius and variable bank angle as in the algorithm of 
ref. (2) . Constant radius turns simplify the computational 
problem. The present study Improves the fuel efficiency of 
the algorithm in ref. (2) for certain types of capture 
trajectories by generating the finei turn of the capture 
trajectory at near constant bank angle. The non-circular 
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turns are computed by a real-time algorithm which uses the 
algorithm of ref .(2)ln an iterative loop. In addition, 
significant chavacterlstics of the optimum airspeed profile 
during turns, determined in ref. (5), are incorporated in the 
synthesis of the airspeed profile. These refinements increase 
the computational complexity of the capture algorithm to a 
modest degree. 

In addition, a more sophisticated algorithm is described, 
which provides a much closer match to the optimal solution for 
long-range flightpaths. Although this algorithm is too time 
consuming for implementation in an onboard computer and is not 
capable of handling short-range flightpaths, it is thought to 
represent the limit to which an on-boar 5 algorithm may be 
improved for the long-range problem. 

Finally, the fuel consumption used in flying trajectories 
computed by tht*. re i-.ime algorithm and the more sophisticated 
algorithm are comp-^ed with the fuel consumption that would 
occur if the aircraft were to fly the trajectories defined by 
the extremal solutions to the optimal control problem. 


In this initial study the effects of altitude and winds 
on the trajectories were ignored. The impact of these 
assumptions on the algorithm is currently under study. 


REAL-TIME ALGORITHM FOR THE CONSTRUCTION OF FUEL CONSERVATIVE 

CAPTURE ’'LIGHIPATHS 

Overview 


Two major steps are involved in deriving the real-time 
algorithm. First it is necessary to identify and briefly 
review the salient features of certain extremal fuel optimal 
trajectories determined in ref. (5). Then it is necessary to 
approximate these characteristics for use in the algorithm. 
Since relatively long approach flightpaths (over 10 miles) are 
the rule, rather than the exception, the characteristics of 
long flightpaths generated, as in ref. (5), were examined. 
Such paths contain a beginning and a final turn, connected by 
a line segment of extremely small curvature, which is almost a 
straight-line. The speed profile for such a path is as 
follows; To achieve a tighter initial turn there is often a 
minor speed reduction followed by an acceleration, which 
continues in the almost straight-line segment. In the 
straight-line segment the speed of minimum fuel per unit 
distance is approached but never achieved. Deceleration at 
zero thrust toward the landing speed usually takes place in 
the final turn, a considerable portion of which is flown at 
the maximum bank angle limiv. If the final turn is large the 
deceleration is often followed by a brief acceleration. 
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It is shown in ref. (5) that the controls, the bank angle 
and the thrust, are continuous in time. To generate and store 
buch control commands would take a large amount of computer 
time and storage. Therefore, in the real-time algorithm 
developed here, we shall be satisfied with step commands in 
the controls. Throttle control will consist of maximum thrust 
in the acceleration segments of flight, thrust to equal drag 
in the constant speed segn.ents, and zero tJirust in the 
deceleration segments. Also, we shall be satisfied in this 
algorithm with flightpaths that are constructed of a 
combination of straight-lines and sequential circular arcs, 
which approximate constant bank-angle turns. For simplicity, 
acceleration segments will occur only on the straight-line 
segment of the path. That means that the initial turn is a 
constant radius, constant speed turn. The bank angle at the 
start of each segment will be either zero or 30°' which is the 
maximum allowable bank angle for a commercial jet liner. This 
bank angle gives the tightest possible turn for the given 
speed; indeed in the optimal trajectories (ref. (5)) tfie turns 
quite frequently do occur at the bank-angle limit. The speed 
profile and the shape of the flightpath are interconnected, 
since the turn radii depend on the speed of the turns. 

Flightpath Development 

Turning Segments . - To construct a flyable flightpath which 
includes a sequence of circular arcs to approximate a constant 
bank-angle turn, the equations for the circular arc segments 
of the flightpath are first developed. The use of these 
segments to approximate the constant bank- angle turn is then 
discussed. 

During constant radius decelerating turns the bank angle, 

♦ , is related to the instantaneous velocity by 

R ■ v^/(g tan^} (1) 


where R is the radius of turn. R is chosen so that, 
maximum bank angle of 30° occurs at the maximum speed 
From F.q.(i),(7) and (8) of Appendix I with T » 0 we have 


dv/dt - -{kjV^ t k^/v^) 


(2) 

where k^ and 

k^ are 




kj • g(kj ♦ k2/(gV)/W 


(3) 


k^ ■ 9 ^2/W 


(4) 


Eq (2) Includes the effect of the varying bank angle while 
decelerating on a constant radius turn. 
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To calculate the end velocity v, from the start velocity v, on 
a circular arc of angle 0 or vice versa* dt inKq.(2)is replaced 
by its equivalent ds/v. Integration gives 


((V, * 


(5) 


where v„ 
m 


is the minimum drag velocity 


m 


(k4/k3) 


1/4 


( 6 ) 


In specific applications the aircraft model may be more 
complex, in which case numerical integration of the equations 
will be called for as in the algorithm of ref. (2) . 

When the final speed and the angle of the circular arc, 

0, are given, and the radius of turn and initial speed are the 
unknowns (ref. (5)) has to be solved by successive 
approximations, since the unknown v^ cannot be isolated. 

The constant bank-angle turn is approximated by a 
sequence of circular arcs of angle 0 with varying radii. If 0 
is too large the aircraft will fly a decelerating turn that is 
less tight than it theoretically could, and the resulting 
riirl'tpath would normally be less fuel efficient. In 
addition, upon transfer from one circular arc to the next a 
large transient in bank angle would occur. On the other hand, 
if we choose the angle 0, too small, computation time and 
computer storage increases in the airborne computer. As a 
compromise a 0 = 30*^ segment was chosen. Usually, the final 
speed and heading were chosen to be 180 knots and 0 degrees 
respectively. Starting from this end condition, by using the 
relationships of equations (1) and (5) , the decelerating 
spiral turn can be constructed as shown in figure 1. 

The Horizontal Trajectory . - The construction of the 
horizontal trajectory is carried out with the help of the 
horizontal capture trajectory algorithm described in rof. (2) 
and (6) . This algorithm was originally designed to provide 
minimum time, (and near minimum fuel) , constant speed 
trajectories. It will construct a trajectory from a given 
initial position and heading to a final position and heading 
by means of a turn along a circular arc followed by a straight 
line (or, less frequently, by a second circular arc) and a 
subsequent circular turn. Depending on the end conditions, 
the algorithm delivers between a minimum of two to a maximum 
of six different paths, of which the shortest path is chosen 
unless a specific type of flightpath is specified (e.g. the 
final turn must be a right turn) . As input, the algorithm 
must be given the starting and the final turn radius, and the 
user must make sure that the desired speed profile can be 
accommodated by this flightpath. As shall be shown, this 
algorithm can be used sequentially to generate the approximate 
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Figure 1. - Construction of a Decelerating Turn Flight Path were waypoint h 2 

defined by position (x^^, y^), heading if , nominal speed at the 
wa:*po1nt, radius of turn nominal control settings ■ 
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constant bank -angle deceleration dlscusaec' in the previous 
section. 

For simplicity, the Initial turn radius is kept constant. 
Assume that the starting and final position, heading, and 
speed are given as sho\ n in figure 2a. For the moment also 
assume that deceleration is required in the entire last turn. 
From Eq. (5) we solve for v , the initial velocity of the nth 
segment where v « v-, and using v, and v ,, we solve for 
and R from Eq. (1) . Application of tnia algorithm (ref. 
(6)) ^irerates the path of figure 2b. This path could not be 
flown, however, since the turn from waypoint n - 1 to waypoint 
3 is tighter than the speed permits. So we use the conditions 
at point (n-1) as the final position speed and heading and 
repeat the above process. This results in the flightpatn of 
figure 2c. We repeat the process until the final circular 
segment is 30°, as shown in figure 3d. The resulting speed 
profile is shown as solid line in figure 2e where the 
deceleration starts in the straight-line segment. Additional 
fuel savings can be had if changes in speed are allowed such 
as is shown in a dashed line in figure 2e. This is further 
discussed in the next section. 

Occasionally, for short flightpaths, application of the 
algorithm in ref. (6) to approximate a conr.tant bank -angle 
turn in the manner described above may not result in the 
circular segments all going in the same direction. Instead 
the algorithm may detemdne that a turn to point n-i in the 
opposite direction is raore optimal. Because of this, the 
real-time algorithm developed here was modified to compute 
both the path and the fuel that would be used if the final 
turn were to be constrained to a single direction and to 
compute the path and the fuel that would be used without this 
constraint. The path using the lt.ast fuel is then selected. 

In general, the final deceleration in the flightpath may 
begin in the initial turn, the straight-line segment or the 
final turn. The algorithm, which is described in greater 
detail in the appendix, provides for these different cases. 

If the flightpath is too short to acconmodate the required 
speed change, the flightpath is rejected and the next longer 
flightpath is chosen from the available set. 

Speed Profile for the Straight-Line Segment . - The flight 
path developed in the last section has a straight constant 
speed segment. If the speed is allowed to exceed the initial 
speed v. (e.g., 250 knots for the terminal area) further fuel 
reduction may be obtained by varying the speed in the 
straight- line segment. For a design criterion we turn to the 
minimum fuel extremals for straight-line flight, an example of 
which is shown as solid line in figure 3. The constant speed 
which is approached but never reached is the speed for minimum 
fuel per unit distance v *, and the thrust is continuous but 
is zero in the latter pSrt of the deceleration segment. 
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Figure 2. - Construction of a Fuel Conservative Capture Flightpath . 
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Figure 3. - Optima, and Suboptimal Speed and Thrust Profiles for a Straight 


Line fllshtpath. 


Equations (2) and (5) apply to decaleration on a straloht-llna 
path by setting the second term In Eq. (3) to zero since R - 
O’. We approximate the optimal profiles by the dashed line of 
figure 3. This approximation Is, In fact, the optimal 
solution when c, ■ 0 In Eq. (9) of Appendix I. It Is 
Identical to theT one used In ref. (2) . Such a solution Is 
simple to Implement, and, as will be shown In the results 
section, comes close to the true optimal solution. 


A MORE SOPHISTCATED ALGORITHM FOR LONG RANGE FLIGHTPATHS. 

The previously described fllghtpaths consist of only 
straight-line and circular elements. Therefore, they are 
completely definable by a series of waypoints. Also, an 
airborne digital guidance system exists at the Ames Research 
Center which provides guidance signals along such paths. 
Hence the guidance problems for such paths have already been 
solved completely and selected fllghtpaths could be flight- 
tested quickly if desired. In contrast, the turning segments 
of a flightpath with true constant bank-angle turns cannot be 
represented in such a manner. New equations for a cockpit 
path display to the pilot and path error nulling guidance 
equations would have to be developed. Still it is of interest 
to determine how much improvement a constant bank-angle turn 
can provide w^en coupled with a speed profile that more 
closely simulates that of the optimal solution. It is to be 
noted that, to keep programming simple, this algorithm was 
intended for long fllghtpaths only, (with a straight segment 
of at least 10 nautical miles) . 

The speed profile chosen is that alven In figure 4. If 
the initial turn is greater than lOO" then we decelerate 
initially without going below v . . Until the turn remaining 
is 100°. We thM iccelerate to«^d v *. If the initial turn 
is less than 100° we immediately accelerate. The final turn, 
which is generated in the bacJcward direction, the angle of the 
decelerating portion of the turn, , is a function of the 
final speed according to 

6.732 - 0.01179 (7) 

Where ^ is in radians and in ft/sec. 

This expression was obtained as an approximaticm of certain 
extremal trajectory characteristics. This makes 
degrees for a final speed, v., of 250 knots, and * 190 
degrees for a final speed of 180 knots. The remaining turn is 
flown at constant speed cx maxim u m acceleration so that v. is 
met at the end point. For final turns less than tXm ab<^e 
numbers there occurs deceleration only. 
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Having defined the speed profile along the path, we turn to 
the path construction. Again the capture algorithm of ref. 
(6) is called to determine turn directions and approximate 
subtended emgles of the initial and final turns and 'll- 
of figure 4. It must be remembered that the algorithm of rer. 
(6) delivers circular turns, while the desired turns are 
constant bank- angle turns. We select the larger of the 
turns 'J' , or 'i'- to calculate the path of a true constant 
bank-angle turn by numerical integration of the equations of 
motion either in the forward or backward direction while 
following the speed profile. Since the 4'*s from the algorithm 
in ref. (6) are only approximate, a path corresponding to only 
one-half of the indicated turn angle is calculated. This 
process results in a new end point for the remaining path to 
be developed. Then the algorithm of ref. (6) is called again 
from the new end point. The above process is repeated until 
the remaining initial and final turns are less than one 
degree. At this time the less than one degree circular turns 
determined by the algorithm in ref. (6) are flown. The 
development of the speed profile for the remaining straight- 
line path is identical with that of the algorithm. In 
addition to the technical difficulties of this algorithm 
described in the introductory paragraph we note that it is 
also computationally lengthier. It is therefore not suggested 
for onboard mechanization but is presented here as a limit 
that may be approached by the algorithm described earlier if 
one would include the more complex spe^ profile and r^ke the 
circular turn segment smaller, (e.g. 15® instead of 30°) . 


EVALUATION PROCEDURE 


The real-time algorithm always constructs a flyable 
flightpath given any set of initial and end conditions. The 
more sophisticated algorithm constructs a path provided the 
end conditions are separated by 10 miles or mors. In 
contrast, when extremals are generated from the optimal 
control solution described in Kreindler, and Neuman (ref. 5) 
only the end or initial condition can be specified in advance, 
not both. Therefore many extremals had to be generated to 
obtain a set of flightpaths which were judged to be 
reasonable test cases for terminal area approaches. Further, 
the exfemals were divided into two classes. Relatively long, 
gently turning flightpaths such as would normally be expected 
on approaches to an airport, and short tightly turning flight- 
paths, that may represent emergency maneuvers. The fuel used 
on the extremal was compared with the fuel used on the 
equivalent path generated by the two algorithms. These 
comparisons are not intended as a complete statistical 
analysis, which would be very complex and is not warranted at 
this time. 

At present, there is no general method to prove the 
(global) optimality of a given extremal. Therefore, some of 
the extremals that were picked for comparison were shown to be 
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not globally optimal by constructing a more fuel efficient 
path via the algorithm in conjunction with a "tree search". A 
more complete discussion of this problem is covered in 
ref. (5) . 


RESULTS AND DISCUSSION 

Evaluation of the Algorithm . - Figure 5 shows the initial 
statei for S6 relatively long approach paths that resulted 
from the computation of extremals (see ref. (5)). Two flight- 
paths are drawn in for illustration. Also shown are 
comparisons of the fuel used when flying an extremal 
trajectory with the fuel used when flying the trajectory 
generated by the real-time algorithm and the more 
spohisticated algorithm. The real-time algorithm provided 
paths that used on the averuye only 1.53 percent or 10.2 lbs 
more fuel than the optimal solution. The more sophisticated 
algorithm further decreased this to an average of 0.57 percent 
or 3.7 lbs. For long, speed-constrained flightpaths (v „ = 
250 knots) (not shown) the suboptimal paths perform even 
better in comparison with the extremals, because the extremals 
have a long, almost straight-line segment with constant speed 
v(t) = V ., which, of course, is easy to approximate. 

To get an idea of the effectiveness of the real-time 
algorithm during the straight-line segment, we study long 
straight paths with an initial speed of 250 knots and a final 
speed of 180 knots. Comparing two possible suboptimal speed 
profiles a flight with a constant speed of 250 knots followed 
by a final deceleration to 180 knots, with a flight that 
initially accelerates toward v * (the optimum cruise speed) 
we find that the latter is 0.9 percent (or 16 lbs) better for 
an 8 nautical mile path and 5.18 percent (or 21.8 lbs) better 
for a 16 nautical mile path. Hence, the initial acceleration 
pays off if the path is sufficiently long. When the sub- 
optimal paths which contain the initial acceleration are 
compared with the extremals we find that they are 1.19 percent 
(or 5.61 ]bs) off the optimal for a 12 nautical mile path and 
only 0.18 percent (or 3.4 lbs) off the optimal for a 71 
nautical mile flightpath. This agrees with our general 
observation; the longer the flightpath, the closer the real 
time algorithm solution will compare with the optim.al flight- 
path. 


Figure 6 shows the end states and fuel consumption for 33 
short approach paths. Three extremals are drawn as dashed 
lines for illustration. The Initial states were obtained by 
backward integration of extremals from a final state. The 
data of figure 6 are ordered with respect to the initial 
speeds, and again, within each category, with respect to the 
fuel used on the flightpath. The flightpaths that are marked 
as coasting trajectories had zero thrust for the entire 
flight. Note that some fuel is consumed in these cases, 
because the engines are not shut off even though thrust is 
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Figure 5. - Initial Conditions and Test Results for 28 Long Approaches. 
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retarded to zero. (This is reflected by c in Eq. (9) of 
Appendix I) . The others had speed reductions below the final 
speed and a buildup of thrust towards the -.nd of the flight. 

Figure 7 shows a horizontal track of some of the extremals 
nuirbered to correspond to the path identification number of 
column 1 in figure 6. Figure 7 is merely intended to give an 
appreciation of the complex patterns of turns of these 
extremals. Many of the highly turning extrem.als can be 
thought of as path streching maneuvers, to dissipate the 
kinetic energy while decelerating to the landing speed. It is 
clear that a simple algorithm that provides for two, or, at 
most, three turns cannot come close to generating flightpaths 
v;ith similar fuel consumiption. 

The real-time algorithm often delivered paths that were 
shorter than the extremal, but they were too short to complete 
the deceleration that is required. Hence the next longer 
flightpath developed by the real-timie algorithm was chcsen, 
v;hich then required substantially more fuel (see the column 
labeled "Real-Time Algorithm" in figure 6) . VJhile the real- 
tine algorithm, shows a fuel consumption up to several tim.es 
that of the extremal, the algorithm always delivers a 
realizable flightpath. Also, we note that in two cases the 
algorithm, delivered a path that was more fuel efficient than 
the extremal, ,This brings us to the necessity of an attempt 
to check ex^i^eftals for global optimality before using them, as 
test cases for a suboptimal high speed algorithm. This was 
done by a "tree search" which found seven additional 
ncn-cptimal extremals. As previously mientioned, this probleru 
is discussed in ref. 5. 

SUMMARY AND CONCLUSIONS 

The real-tim>e algorithm always delivers a realizable flight- 
path. It is fuel conservative and near fuel optim.um for most 
flightpaths occuring in the term.inal area. Hov?ever for the 
few short em.ergency maneuvers, which primarily consist of 
highly turning paths flown with zero thrust the algorithm v;ill 
deliver a flightpath that will be far from fuel optim.um. 

The more sophisticated algorithm presents an upper limdt to 
the performance for the onboard algorithm during long range 
flightpaths. This is, for long flightpaths, the algorithm miay 
be im.proved from being within 1.5% of the optimal solution to 
C.6% of the optimal solution. 

The variable turn radius feature of the real-timie algorithm 
saves fuel but adds computation time when compared to a 
previously developed algorithm. Before this algorithm can be 
incorporated in an onboard guidance systemi it must be expanded 
to include the effects of altitude, winds and navigation 
errors . 
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Figure 



APPENDIX I 


AUOAPT MXEL 


Ihe aircreift model used is an e^sproodmation of that of a Boeing 727 
aircraft. Since it is identical with that in ref. (5). it is restated 
here without detailed developnent. The bank angle is limited to 


- 30° < <30° 

The thrust T is limited to 
0 T <_ 30000 lbs 

For horizontal fli^t vertical lift equals weight (fig. 8a) 
L cos = W 

The point mass equations of motions cire as folla*;s: 

The ground velocity catponents aie (fig. ) 

X = V cos ft/sec 

y = V sin ^ ft/sec 
The heading rate is 

X = -g tan (|) rad/sec 

V 

The acceleration along the velocity v-ector is 
V = g (T - D/TV ft/sec 
The drag is 


D = v + 


^2 (1 + tan^ 4i) lbs 


The fuel flow is 


V 


2 

f = c ->• c, T + c» T"^ lbs/ sec 
o 1 2 


( 1 ) 


( 2 ) 


(3) 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 


The constants for a typical ccrmercial jet aircraft in the velocity 
range of 150 - 350 knots are; 


k^ = 0.02808 lbs sec^ ft“^ 

k 2 = 606,055,000. lbs ft^ sec"^ 

c = 0.80833 lbs sec”^ 
c 

Cj^ = 0.000150694 sec 

= 5.4 X 10 sec ^ lbs ^ 


g * 32.2 ft/sec 
W =150000. lbs 


-1 
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L cos 0 



a) REAR VIEW OF AIRCRAFT 



» X 


b) TOP VIEW OF AIRCRAFT 


Figure 8. - Problem Coordinate System. 
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APPENDIX II 


SUBROUTINE DESCRIPTICN FOR IHE FEAL-TD® ALOORITWl 


Subroi’tijie Linkage - 


Ceillirg Program 


MAIN 
I£\;EL 3 
FDEX3L 


SUBROUTINES CAU£D 



SLFP 


X 
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MAXDEC 



SUBROUTINE I£VEL 3 


Given: VS ® starting velocity, V4 ■ final, velocity, starting position and 
heading, and £in£d position and heading. 


Desired: Oonstnx± a flyable trajectory in essentially real tine, for ary 
two end conditions (VS ^ V4) \Mch is a reasonable approxination to a 
nlniiTum fuel fli^t path. 


Such approxinate paths \asualiy have an initial turn, possibly with variable 
radius, a straight-line flight segment, and a final variable radius turn. 
For fuel efficiency, the required deceleration is perfomed as late on the 
flightpath as possible, and turns are made at the maxirun bank angle 
allowable for the given speed. Depending on the length of the path and the 
speed differential between starting and final velocity, deceleration may 
occur in all three fli^t segments. The final variable radius turn is 
approximated ly a series of constant radius turns, each with a heading 
change of 30 degrees or less. 


1. Calculate the turn radius Rl of the initial turn from VS the starting 
speed. Check if final speed is greater than initial speed. If yes go 
to 25. 


2. Call VSTART. Calculate the turn radius R2 and the initial speed VI for 
a 30 degree circular arc porticn of the final turn from the final 
speed V2, assuning maximun deceleration (thrust 0) and iraxijiun bank 
angle at the beginning of the segment. 


3. If VI calculated in step 2 is greater than the initial speed VS, because 
the actual remaining po^on of a circular arc is smaller than 30 
degrees or because deceleration begins in the final turn after an 
initial constant speed turn, the turn radius calculated would be larger 
than required and will be set to the maximum vedue R2 bR 1, (Rl determined 
in step 1) . 


4. Call CAFT(6) vhich constructs the shortest flightpath using Rl and R2 
as initial and final turn radii iron the specified initial to the 
^lecified final point and heading, using only circular turns and a 
straight line. (If desired CAFT can be made to consider only a final 
right or a final left turn) . 
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5. OiecJc if final turn is less than 30 degrees (turn oonpleted) . If yes go 
to step 10. 


6. Check if this 30 degree segnmt is VI greater than VS, that is could we 
decelerate from VS to V2 in less than a 30 degree turn. If yes go to 
19. 


7. Call MAXD. Frcn the xadius R2, and the initial and final ^^eeds of the 
30 degree turn segment VI and V2 calculate the time to fly, fuel, and 
distance. 


8. Call NEVUS’. Calculate the position of the beginning of the turn 
segment, which defines a waypoint anA forms the new enc^int for the 
continu^ backward development of the flight path 


9. Go back to step 2. 


10. Call MAXD. First segment of the final turn. Using the remaining turn 
angle and R2 and V2 calculated earlier, find the speed at the beginning 
of the fir.al turn and the length of the first circul 2 u: azc. 


11. If VI calculated in step 10 is larger than the naximtm speed VS the 
deceleration begins in the first segment of the final turn, go to 19. 

12. Call itie deceleration starts before the last turn. Depending on 

whether the trajectory found last by CAPT(6) was a tum-straight-tum or 
a tum-tum-tum cme the distance S is calculated that it vrould take on 
the iriddle segment of flightpath to decelerate frem the inaxiimir 
velocity VS to the speed at the beginning of the final turn. 

13. If the length of the middle segnent, D, calculated the last tine step 4 
was executed is longer than S fron step 12 the acceleration begins in 
the irlddlc segment, else go to step I.T. 


14. From the ratio SA> and the location of the end^ints of the middle 
se^nent the position of the beginning of deceleration is calculated. 
(This calculation is different for a strai^t>li^ and a tuming~oenter 
secpient) . Go to 2!^ . 
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15. ihe deceleration begins in the first turn. Calculate the qwed at the 
end of the first turn, VI, (vhich is also the beginning of the middle 
segment) and the fuel and time for the middle se^nent, which is 
oonpletely flown at zero thrust. 


16. Chedc if the first turn is long enough for the required deceleration, by 
calculating the distance SD it would take to decelerate from VS to VI 
with a turn radius R1 and oonparing it to the actually available 
distance c£d.culated in step 4. If okay go to 18. 


17. Path required is longer than available. Print that path stretching is 
required. END. 


]8. Path is long enough to have a constant speed segment. Calculate vhere 
constant speed segment begins and to to 23. 


19. Deceleration begt^ in the last turn. Calculate the time and fuel for 
the deceleration (MAXD) and position (NE5CWP) of the beginning of the 
deceleration tui.;. 


2C. Calculate the remaining constant speed flight path CAPT(6) and calculate 
fuel and time on the initial part of the last turn. 


21. Calculate fuel and time for the middle segment (straight or turning) . 


22. Check if middle segnent is strai^t or turning. If strai^t go to 24. 


23. Call FUELC. Calculate fuel and time for the initied constant speed 
turn. EUD. 


24. Ihe center segment is a strai^t line. Calculate how nich fuel can be 
saved by accelerating toward the optinun cruising velocity at the 
beginning of the segment and decelerating back to the nominal velocity 
at the end. If this maneuver is permitted (no speed limiting) then the 
possible savings is subtracted from the total fu^ for the fHghtpath 
under oonaideration. Go to 23. 


25. Vhen VS1V4 first and last turn are flown at at constant speed with 
mininun radii for the given speeds. Fly the straight line with initial 
acoeleraton towards the speed for maxiirun distance per unit of fuel 
final deceleration towards V4 (call SUP) E2C. 
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FDECEL (XI, Yl, Hi, VS, V4, J, L) 


Qilculates for an initial ri^t (L ■ 1) or l«ft (L ■ -1) turn tuooaMive 
ooordinates for a 0 thxuat deoelerating tum^ including naw zadiua and 
new starting sp eed, total tine and fuel used, and prints out the data, 
vdhere J is the nuiter o£ 30 degree incronents. 


Purp^ t It) peonit an initial deoelerating turn to the left at right before 
developing the remainder of the fiightpath in the l£ViI, 3 svtiroutine. 

''1 



Given the end velocity V2 and the turn angle 6 the subroutine calculates the 
starting veloci^ VI and the radius of turn R given that the initial bank 
angle is 4 *= 30~and a power off decelerated turn is perfonned. 


MAXDBC 



VS ■ Starting Vel, 

VF - Final Vel. 

D - Distance 

R ■ Radius ctf turn (P » 0 for 
strai^t line). 


Givoi two of the three quantities VS, VF, D and R the subroutine calculates 
the third quantity via Numerical integraton for a power off deceleration of 
the aircraft. In addition it calculates the time and^fuel for the.naneuver. 
MAXD same as MJOTOC but uses analytical expressions for Cic calculations. 


XlVlh, 



Calculates heading ai d iosition (X, Y) of the next viaypoint from the 
position and heading or the old %Mypoint as a function of direction of turn, 
turn angle 6 , direction of evaluation (get X 2 Y 2 
frem Xj, • 


PPOT’ 

This subroutine converts waypoint variables to desired dimensions for 
printout, prints out the variables and transfers data to calculate ncsxt 
waypoint ooordinates. 
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OffTCMX) 

OCMCN/RjEVSTM/r(20), JT 


The subroutine capt will generate from a given position (x. , y., i>.) to a 
final position and heading (x., y., il'.) a>id initial and fij^ Cum radius P. 
and R. the mininun length fllghtplth from the six possible candidates ahe^n 
on the next page. Fcac the purpose of finding an ipproodmation of a mininun 
fuel path vari^ types of can^dates may be exclu^ from the search 
the nGbinun length path. 

^K « 4 exclude the Turn-turn-turn trajectories from the search (fig. 9.5 
and 9.6 next page) 

MK ■ 6 include all types of trajectories 
From the types of trajectories not excluded by MX in addition exclude 


JT - 1 


none others 


JT = 2 all that have a left turn last (RSL, LSL, l^) 

JT 1, JT 2 all that have a right turn last (PSR, I5R, RUl) 


Types of fi lehr paths considered by the algorithm for generation of 
torlzcntal capture trajectories from a given position and heading to a final 
position and heading. At least two and at most 6 solutions exist when R. 
and R^ are given. The shortest path from the exiucing solutions is chos^. 
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SLFP 


This subroutine calculates time and fuel It ta)ce8 for a straight 'line flight 
path of distance D, starting speed VS and final speed VF. The strategy that 
is used would be fuel cptijnum if fuel flow rate were linearly related to 
thrust (f * k T) . 

Strategy: Accelerate with max thrust towards the speed at v^ch min. 
Whether it is reached or not, start decelerating with 0 thrust at the proper 
point in time so that the firal speed is reached at the distance D. 



(B) 



This is acccmnlishecl as follows. A distancje D, to accelerate or decelerate 
frcm V to is calculated and cxnparod with D , the available distance. 

If D if smaller than D no flightpath exists toHmeet the ocnditicns. If 
D if larger than D, the ranainlng path is generated via nuterical 
integration of the accelerating and decelerating segments simultaneously 
until either the loitering speed is reached (A) or until the path length D 
has been achieved before reaching the loitering speed (B) . In c:ase of A a 
constant speed segment is added to achieve the carplete distance D^, and in 
case of B the calc:ulations are finished. 
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